o 
o 

(N 



> 



Perturbation Theory for Population Dynamics 



Francisco M. Fernandez 



o ; INIFTA (Conicet, UNLP), Division Qmmica Teorica, 

Q ; Diag. 113 y 64 (S/N), Sucursal 4, Casilla de Correo 16, 

O ■ 1900 La Plata, Argentina 

^ ■ E-mail: fernande@quimica.unlp.edu. ar 

February 2, 2008 

Abstract 



We prove that a recently proposed homotopy perturbation method 
for the treatment of population dynamics is just the Taylor expansion 
of the population variables about initial time. Our results show that 



1 Introduction 



, this perturbation method fails to provide the global features of the 

I ecosystem dynamics. 

(N 

o 

Recently, Chowdhury et al [1] proposed the application of homotopy-perturbation 
I method (HPM) to simple models of population dynamics and obtained ap- 

^ ■ proximate solutions in the form of perturbation series. The authors show 

that their approximate analytical results agree with the numerical solution 
of the problem and appear to suggest that the HPM series are convergent. 

A straightforward inspection of such series reveals that they are merely 
Taylor expansions of the time variable. One does not expect such a local 
approximation to provide a reasonable description of the dynamics of non- 
linear systems, except in the neighbourhood of the initial conditions. Singular 
points appear spontaneously in nonlinear systems and move around the com- 
plex plane as the initial conditions vary [2] which makes unlikely that the 
time series are valid for all time. 
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The purpose of this letter is to investigate the range of utihty of the ho- 
motopy time series to provide useful insight on the dynamics of population 
models. In Section [2] we show that the HPM of Chowdhury et al |lj always 
leads to a Taylor expansions of the solution of the nonlinear system about ini- 
tial time. In Section [3] we analyze the exactly solvable one-dimensional prob- 
lem already considered by Chowdhury et al [Ij. In Section H] we study their 
population model for two species [T] and an exactly solvable two-dimensional 
dynamical model. Finally, in Section [5] we summarize our results and draw 
some conclusions. 



2 Perturbation method 

Population models give rise to differential equations of the form 

x(t)=f(x(t)), x(0) = Xo (1) 

where x is a vector of the n variables xi,X2, . . . ,x„ and f(x) is a vector- 
valued function with components /i(x), /2(x), . . . , /„(x). We assume that 
this vector-valued function is continuously differentiable. 

The HPM proposed by Chowdhury et al pQ is equivalent to a straightfor- 
ward perturbation theory based on the modified equation x(A, t) = Af (x(A, t)), 
where A is a dummy perturbation parameter. Notice that x(l,t) = x(t) so 
that we set A equal to unity at the end of the calculation. The HPM proposes 
a solution in the form of a series 

oo 

x(A,t) = 5:x(^)(t)A^- (2) 
i=o 

where A is finally set equal to unity as said above. Notice that since x'^^^ (t) = 
the resulting unperturbed or reference model x^^\t) = xq seems to be 
quite poor at first sight. Our results clearly show that it is actually the case. 
The initial conditions for the perturbation corrections are x*^-^^ (0) = for all 
j>0. 

If we define the new time variable r = At, then rfx/c/r = f(x) from which 
we conclude that x(A, t) = x(At) and that this particular implementation of 
the HPM becomes the straightforward time series 

oo 

x(t) = Ex,t^- (3) 

i=o 
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in agreement with the particular results derived by Chowdhury et al [T] . For 
example, the coefficient of ffist order is simply xi = f(xo). 



3 One— dimensional model 

We ffist consider the simple one-dimensional model [1] 

X = x{h + ax) (4) 

where 6 > and a < 0. In this case we have an unstable node at x = and 
a stable one at x = —h/a [2j. The exact solution is 



b>0 



As mentioned above, the HPM series agree with the Taylor expansion of 
the exact solution about t = 0: 

t^xo(2axo + b){axo + h) 



x{t) = xq + txo{axo + h) + 



2 

t^xoiQa^XQ + QahxQ + 6^)(axo + h) 



6 

t^Xo{2axQ + h){l2a^x1 + 12abxo + 6^)(axo + b) 
+ ^ +••• ( ) 

Equation clearly sows that the solution has a pole at = ln[l + 
b/{axo)] and, therefore, the HPM series does not converge for t > \tc\- When 
b = 1, a = —3, and xq = 0.1, then \tc\ = 3.253846656 is much larger than 
the largest time value chosen by Chowdhury et al [Ij, which explains why 
those authors obtained such good results. If, for example, xq = 1, then 
tc = — ln(3/2) and the HPM series is unsuitable for t > ln(3/2) as shown in 
Table [H Notice that these results reflect the fact that the usefulness of the 
HPM depends on the initial conditions. Besides, the HPM series does not 
take into consideration the stable node at x = —b/a and is therefore unable 
to reveal the main features of the dynamical behaviour of the system. 



4 Two— dimensional Models 

Chowdhury et al \i\ also discussed the simple two-species model 

X = x{bi + aux + auy) 
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y = yip2 + aaix + a22y) 



(7) 



and obtained accurate results for the model parameters hi = 0.1, an = 
-0.0014, ai2 = -0.0012, 62 = 0.08, 021 = -0.0009, 022 = -0.001, and initial 
conditions xq = 4, ?/o = 10. However, the time interval considered by the 
authors is too small to give any indication of the evolution of this model 
ecosystem. 

By inspection of the first terms of the HPM expansions 

x{t) = Xo + txo{auXo + ai2?/o + h) + t'^Xo[2alixl + 3aiiXo(ai2?/o + h) + aj^yo 

+ai2yo{a2iXo + 0221/0 + 26i + 62) + bl]/2 + ... 
yif) = yo + tyo{a2iXo + a22yo + b2) + t^yo[aiia2ixl + ai2a2iXoyo + ali^o 

+a2iXo{3a22yo+bl + 262) + 2a222/o + 3a22&2l/o + bl]/2 + ... (8) 

we appreciate that if the model parameters 6j and aij are sufficiently small 
(as those chosen by Chowdhury et al [IJ), then the time series may give 
accurate results for an apparently large time interval. However, these series 
do not take into consideration the critical points of the model equations and 
therefore they cannot reveal the actual dynamics of the system [21 [3] . 

The nonlinear dynamical system ([7]) exhibits four critical points in phase 
space [2], E]: an unstable node at (x, y) = (0, 0), two saddle points at (x, y) = 
(0,80) and {x,y) = (71.43,0), and a stable node at {x,y) = (12.5,68.75). 

Fig. [1] shows that the population moves in phase space from the initial 
condition to the stable node and that the time series is unable to take into 
account this important dynamical behaviour. Increasing the perturbation 
order from four (the one used by Chowdhury et al pQ) to ten just improves 
the accuracy for small time but worsens it at larger time which suggests that 
the convergence radii of the time series are rather too small. In other words, 
the time series do not allow us to study the important population portrait 
in phase space. 

In order to appreciate a more dramatic failure of the time series consider 
the following system of nonlinear equations [3] 

X = —y + ax{x^ + y'^) 

y = x + ay{x^ + y"^) (9) 

It is unsuitable for population dynamics because it allows negative values 
of x{t) and y{t) but has the great advantage of being exactly solvable. Its 
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solutions are 



xit) 



ro cos(6'o + 1) 



yit) 



ro sm(6'o + 1) 



- 2arlt 



^1 - 2arlt 



arctan — 




(10) 



We have an unstable spiral when a > and a stable one when a < |3j . We 
clearly see the pole at tc = l/{2arl) and realize that the time series will be 
completely useless for t > \tc\. 

Fig. [2] shows the numerical and approximate solutions for this model 
when a = —0.5 and Xq = i/o = 2. We clearly notice that the time series 
fail completely to provide a qualitative description of the spiral point and 
thereby of the global details of the system dynamics. 

It is worth mentioning that multiple-scale perturbation theory [21 E] gives 
the exact answer for this model and therefore appears to be a much more 
reliable perturbation approach for nonlinear dynamics. 

5 Conclusions 

We have shown that: 

• the homotopy perturbation method proposed by Chowdhury et al [1] 
is just the Taylor expansion of the solutions of the nonlinear systems 
about t = 0. 

• the perturbation series, and consequently the HPM, is limited to a 
neighbourhood of the initial time determined by the singular point clos- 
est to the origin of the complex t-plane. The locations of the singular 
points of the nonlinear models shift as the initial conditions vary. 

• the HPM does not give an acceptable qualitative description of the 
most important features of the evolution of the dynamical system in 
phase space. In this sense, the HPM is by far less useful than the 
standard linearization which is also a local approach [21 [3]. 

In principle, other implementations of the homotopy perturbation method 
may be more suitable for nonlinear dynamics. For example, we may choose 
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the linear approximation about the critical or fixed points [21 [3] as unper- 
turbed or reference model for the subsequent application of perturbation 
theory. 

Since homotopy perturbation methods have become quite popular and 
are currently being applied to a wide variety of fields [H H] (and references 
therein), present results become important because they suggest that a more 
careful scrutiny of the approach's performance is required. 
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Table 1: Logarithmic error log \ {Exact — Approximate) / Exact\ for the 
fourth-order time series for model (jl]) 



t 


Logarithmic Error 


0.1 


-3.14 


0.2 


-1.66 


0.3 


-0.798 


0.4 


-0.193 


0.5 


0.273 


0.6 


0.651 


0.7 


0.968 


0.8 


1.24 


0.9 


1.48 


1.0 


1.69 




Figure 1: Numerical (squares), fourth-order time series (solid) and tenth- 
order time series (dashed) curves in the phase plane for model (^^. 
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Figure 2: Numerical (solid) and fifth-order time series (dashed) curves in the 
phase plane for model ([9]). 
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